Quantum simulation of a triatomic chemical reaction with ultracold atoms on a 
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We study the scaling and coordinate transformation to physically simulate quantum three-body 
coUinear chemical reactions of the type A+BC — >■ AB+C by the motion of single ultracold atoms or 
a weakly interacting Bose-Einstein condensate on an L-shaped waveguide. As an example we show 
that the parameters to model the reaction F+HH H+HF with lithium are at reach with current 
technology. This mapping provides also an inverse scattering tool to find an unknown potential, 
and a way to transfer the knowledge on molecular reaction dynamics to design beam splitters for 
cold atoms with control of the channel outcome and vibrational excitation. 
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Introduction. Ultracold atoms and ions are rela- 
tively easy to isolate, prepare, manipulate and detect 
by means of highly controllable operations that preserve 
their quantum coherence in the time scale of processes of 
interest. They have thus become natural candidates for 
performing physical, rather than numerical, simulations 
in which the effective Hamiltonians governing their dy- 
namics can be made equal to the Hamiltonians of very 
different, simulated, quantum systems. These simula- 
tions are thus based on a formal analogy and may pre- 
dict the behavior of the simulated system under condi- 
tions hard to realize and/or calculate in the original one. 
The simulating system may also be interesting in its own, 
beyond the parameters relevant for the simulation, and 
lead to genuinely new phenomena and applications 
This opens exciting perspectives for many-body physics 
0, and also for few-body systems. 

In this letter we show that this approach can be ap- 
plied to molecular dynamics and chemical reactivity by 
studying the analogy between reactive collinear three- 
body chemical reactions and the motion of a single cold 
atom, or possibly a weakly interacting condensate, on 
a potential surface designed by a magnetic or optical 
waveguide. We put the emphasis on the chemical re- 
action, but the same procedure may shed light also on 
non-reactive collisions. What we propose and what is 
facilitated by ultracold atoms is basically a quantum dy- 
namical version of the rolling ball analogy of chemical 
reactions with the ball ensemble substituted by a conden- 
sate or an ultracold-atom wavepacket, and the mechani- 
cal model potential by a magnetic or optical waveguide. 
Quantum effects are important for state to state (rather 
than averaged) results as well as for reactions involving 
a light atom transfer such as hydrogen. 

Most chemical reactions occur with steric require- 
ments, i.e., a preferred direction of attack. The collinear 
configuration for the reaction path corresponds in many 
"abstraction" reactions involving halogen and alkali 
atoms to the lowest potential barrier and to the preferred 
orientation within a narrow cone of acceptance j3| . More- 



over collinear reactions may be induced by orienting cold 
polar molecules with strong electric fields via the second 
order Stark effect Q. They are also a standard work- 
bench for testing new calculational methods, examining 
the range of validity of several approximate theories, and 
exploring parameter variations over a wide range of val- 
ues, difficult to implement with full 3D calculations. Ac- 
curate quantum calculations involve two mathematical 
coordinates and are still time consuming and specially 
troublesome when heavy atoms or high energies are in- 
volved. The results of interest are usually the branching 
ratios among the channels or the distribution of produced 
molecules among the possible vibrational states. 

Simulation Setting. The collinear chemical reaction 
A-I-BC — )■ AB-I-C, corresponds to the collision of an atom 
A and a nonrotating diatomic molecule BC with the three 
atoms aligned. We assume that the Born-Oppenheimer 
approximation holds and separate the fast electronic and 
the slow nuclear motions. In terms of nuclear masses, po- 
sitions and momenta in a laboratory frame the nuclear 
motion is governed by the quantum-mechanical Hamilto- 
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where V is the effective interaction between the three 
nuclei. The first step is the transformation from the 
"chemical reaction" variables to the atomic "simulation 
variables" . A second important task is to show that the 
required parameters for the cold atom experiment are 
available with current technology. 

Mass-weighted coordinate system. Let us introduce 
the center of mass (CM) coordinate Rcm ■= (mAXA + 
rriBXB + rncxc)/M , where M :— rriA + ttlb + mc, and 
the relative coordinates qi := xb — xa, <l2 ■— xc ^ xb- 
The corresponding momentum operators are Pcm ■= 

pA+pB+PC,Pq,:=^^ = ^PcM-PA, 

'cm- The kinetic energy T 
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is not diagonal in terms of them. To diagonalize T, we 
use the mass-weighted coordinates [1] 
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with mass factors a 
yJmc{mB 
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mA)/M, tan/? — ^ msM / (m attic), and 
scahng parameters m and / that we can choose freely. 
The corresponding momentum operators are 



Pc 



h d 



la\ 



i dQi 
h d 



rriB + mc 
~ hi sin j3\fffi 



niB + nic 



niA 



-PA +PB +PC 



1 



1 



-PC 



i dQ2 \mc' TUB 

and the kinetic energy T takes the form 
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The connection between the simulation variables 
{RcAi,QiiQ2} and the chemical reaction variables 
{xa,xb,xc} is given by 
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In the following we ignore the trivial center of mass 
motion and assume that the potential does only de- 
pend on the relative differences between the parti- 
cle positions. Then the corresponding time-dependent 
Schrodinger equation associated with the Hamiltonian 
([T]) in the new variables is 
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where we have set r — t/P and 



* + V^q(Qi,Q2)*, (5) 



VQ{QuQ2)=l^Vq{quq2)^PV[xA,XB,Xc). (6) 

Equation ([5]) is the important result for the simulation, 
and describes 2D quantum motion of a quantum particle 
of mass m on the potential Vq . 

Potential energy surface. We now specify the potential 
surface Vq for the interaction between the three particles 
of the reaction. This might be an ab initio or, more gen- 
erally, a semiempirical potential. Here we assume the 
semiempirical London-Eyring-Polanyi-Sato (LEPS) sur- 
face 
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and qs = qi + q2- Di, Pi and qio are the dissociation en- 
ergy, the Morse parameter and the equilibrium distance 
of the i-th diatomic molecules that we can construct from 
the three atoms. The adjustable parameter A is opti- 
mized for each reaction. In the asymptotic regions, be- 
fore and after the reaction happens, one of the atoms is 
far from the others and the potential energy is the one 
of a diatomic molecule 4]. In the LEPS surface, this is 
given by the Morse function 
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where j = 1 for the products' channel with the diatomic 
molecule AB, or j = 2 for the reactants' channel with 
the diatomic molecule BC. This potential near the equi- 
librium distance qjo can be harmonically approximated 
by 



(9) 



2Dj/3'^ is the force constant. 



where Kj - ^^^^.j 

Applying Eqs. Q to the potential in the asymptotic 
regime where Vq{qi,q2) ~ Vj{qj), we obtain for the sim- 
ulating frame that the energy surface in the asymptotic 
regions of the products' and reactants' channels are, tak- 
ing into account Eq. ©, 
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whereas for the reactants' channel (j = 2) 
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The function xi is a rotation in the (Qi, Q2) plane, so the 
potential Eq. ([T0| is, for the products, simply a rotated 
harmonic oscillator in the (Qi, Q2) plane. In terms of the 
oscillation frequencies of the diatomic molecules Vj, the 
frequencies Vj of the harmonic oscillators in Eq. (|10p are 
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where ^ab and hbc are reduced masses for the diatomic 
molecules. The value of I can be fixed from these last 



equations, so that the potential parameters of the simu- 
lation can be made realistic. 

Initial atomic velocity. To set the initial velocity of 
the cold atom vq-^^ in the reactants' channel we first esti- 
mate the velocities involved in the chemical reaction. If 
the reaction happens at temperature T the rms mean 
velocities for the atom A and diatomic molecule BC 
along a given direction are respectively [ksT /niAY^'^ and 
[ksT / {niB + wc)]^^^, where ks is the Boltzniann con- 
stant. We may then assume 

vc ~ fs, VA - vb = V'aJ^I — ^ + , =y 

V \rriA \/mB + mc J 

(12) 

which, following from Eqs. ([3]), corresponds to the atom 
velocity 




VQ. « 0. (14) 

Example and numerical values. As an explicit exam- 
ple we consider the reaction F-I-H2 — ?> FH-I-H, where 
F-^A, H^B, and H^C, so jua = 3.15 • lO^^^ kg 
and niB = nic = 1.66 • 10"^'' kg. For this particu- 
lar reaction, A ~ 0.164 [1], and the mass factors are 
a = 5.48 • 10"", b = 3.98 • lO'^^ and /3 = 46.45°. For 
the diatomic molecule HF, qio — 0.917A, the dissocia- 
tion energy is Di = 9.609 • 10"^^ J and the Morse pa- 
rameter (3i — 2.242A^^, so the force constant is Ki ~ 
2DiP\ = 966 N/m. Consequently the oscillation fre- 
quency Vl — (i^i//iHi?)i/V27r = 1.246 • 10" Hz whereas 
for HH, (J20 = 0.742, D2 = 7.608-10-19 J, /32 = 1.942l-i 
and K2 = 573.85 N/m so V2 = {K^l mnf^'^ I^t^ = 
1.32 • Hz. 

To simulate the reaction we propose ^Li atoms. One 
advantage of '^Li is that the interatomic repulsive interac- 
tions are extremely tunable with a Feshbach resonance. 
The zero crossing of the s-wave scattering length is the 
shallowest known, so that only modest field stability is 
needed to achieve a non-interacting gas [7| . We thus have 
m = 1.1526 • 10-26^ kg and set I = 6.55 • 10"^. Defining 
the valley depths V, = -DjP (j = 1,2) according to Eq. 
(|6]), the parameters in the asymptotic region for the re- 
actants' channel are = 5.66 kHz, and V2 — 2.4 /iK, 
whereas in the asymptotic region of the products' chan- 
nel, vi = 5.34 kHz, and Vi = 3/iK, The choice of a light 
atom such as lithium is also dictated by the requirement 
of achievable transverse frequencies in the reactants' and 
products' channels with standard techniques (see below). 
To illustrate the scaling of distances and velocities note 
that a displacement of 1 A of the atom F along the reac- 
tants' channel corresponds to a displacement of 7.8 /im 
of the lithium atom according to Eqs. ©. If the reaction 
occurs at room temperature, T = 298 K, Eq. (IT51) sets 
for the lithium atom a velocity vq-^ = 5 mm/s along the 
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FIG. 1: (Color online) (a) Contour map of the potential en- 
ergy surface Q for 1^2 + F HF + H. (b) Contour map 
of the potential for the ^Li atom that simulates the chemical 
reaction. In both cases the energy is in units of the zero point 
energy of the reactants' valley, and the surface is truncated 
well below zero energy (the asymptotic value when all atoms 
are far apart) to better visualize the saddle and reaction path. 



asymptotic region of the reactants' channel. The con- 
trol of matter waves at such low velocities is at reach Q . 
In Fig. [T] we plot the potential energy of the chemical 
reaction H2+F HF4-H given by Eq. dTj to see the 
transformations from the chemical reaction parameters 
{(71,52} into the "laboratory" simulation waveguide on 
which the ^Li atom moves. Note the advanced saddle, 
and the deeper product's valley, responsible for the ex- 
oergicity and the vibrational excitation of the resulting 
HF molecule. 

The experimental realization with ultracold atoms in- 
volves (I) the preparation of a propagating matter wave 
in a guide and (II) the realization of a guide with the 
appropriate shape. A Bose-Einstein condensate, rather 
than the repetition of the experiment with single atoms, 
provides the ideal setting since the fate of the whole quan- 
tum wave packet can be measured in one single exper- 
iment. The propagation of a Bose-Einstein condensate 
into straight magnetic or optical guides has already been 
demonstrated experimentally [l3| • More recently, the 
production of guided atom lasers shows that a large con- 
trol of the matter wave parameters such as the mean 
velocity (5-30 mm/s), the transverse mode occupations, 
the internal state, or the linear atomic density can be 
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achieved pl] - [T^ |. Using different outcoupling mecha- 
nisms, the matter wave can be prepared in the transverse 
ground state In these latter schemes, the dilute- 

ness of the matter wave suppresses the role of interactions 
providing a well-suited system for the quantum scatter- 
ing experiments of interest, without the need of Feshbach 
resonance tuning. The second aspect deals with the po- 
tential modeling to design simple reactive chemical re- 
actions. Different strategies can be envisioned (i) with 
wires sculptured on atom chips by a focused atom beam 
technique 16j , (ii) with adiabatic radio- frequency po- 
tentials |l7l, Il8|. (iii) with high resolution time averaged 
optical potentials "painted" by a tightly focused rapidly 
moving laser beam on a 2D canvas formed by a static 
light sheet [l^. A canvas of 60 /im diameter, and a ra- 
dial condensate thickness of less than 1 as the ones 
realized in , are enough for the spatial range and res- 
olution needed for the simulation, see Fig. [T}d. More- 
over the potential depth can be controlled by velocity or 
intensity modulation, and no decrease in the number of 
condensate atoms is observed after 2.5 s, again more than 
enough for implementing a process of the order of mil- 
liseconds. Reaction probabilities could be detected with 
an in situ high resolution imaging, whereas the coher- 
ent vibrational excitation is measurable after a few ms 
timc-of- flight. A high flexibility in the guide design is 
also provided by combining properly these various tech- 
niques and/or using time-dependent optical or magnetic 
potentials [1^ [2l|. The simplest realization would in- 
volve a crossed red-detuned dipole beams configuration 
in combination with a well positioned repulsive poten- 
tial wall realized by a sheet of blue-detuned laser light 
[23 | . Alternatively, one could study the motion of an ion 
into a well-designed guide. Ultracold ions have already 
been transported in complex structures [H, , but their 
propagation in guides has not been investigated so far. 

Discussion and Outlook. We have worked out the map- 
ping between a quantum-mechanical coUinear triatomic 
chemical reaction and the motion of ultracold atoms on 



a tilted, L-shaped waveguide. As an example we have 
shown that the parameters for simulating the reaction F 
+ H2 — > FH + H using ^Li can be implemented with cur- 
rently available technology. This approach is thus com- 
plementary to other proposals for simulating chemical 



reactions |25| . which are more ab initio and do not need 
any previous calculation of the potential surface or ap- 
plication of the Born-Oppenheimer approximation, but 
require a quantum computation with hundreds of co- 
herently manipulated qubits. This is currently out of 
reach for a reaction like the one discussed. The present 
approach is less fundamental, since it assumes a poten- 
tial surface and the Born-Oppenheimer approximation to 
hold, but also easier to implement. As an inverse scat- 
tering tool, the capability to manipulate the potential 
parameters may be used to fit experimental results of 
the chemical reaction and find the right potential. 

By a straightforward generalization, we could also sim- 
ulate coUinear four-atom reactions by an ultra-cold atom 
in a three-dimensional potential. As a further application 
of the mapping, the vast knowledge and experience accu- 
mulated on chemical reaction dynamics, in particular for 
triatomic systems in the coUinear configuration, is now 
ready to be transferred to design crossed laser beams or 
waveguide bends with different properties. They could 
be used, for example, as control devices for asymmetri- 
cal beam splitting into the channels or for controlling the 
transverse vibrational excitation. An example of this is 
the recent design of an atom diode or one-way barrier 
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